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We perform a calculation of the Sivers function in a spectator model of the nucleon, with scalar and 
axial-vector diquarks. We make use of gluon rescattering to produce the nontrivial phases necessary 
to generate the Sivers function. The inclusion of axial- vector diquarks enables us to obtain a nonzero 
Sivers function for down quarks. Using the results of our model, we discuss the phenomenology of 
transvere single spin asymmetries in 7r+, ir~ , and tt'' production, which are currently analysed by 
the IfERMES and COMPASS collaborations. We find that the inclusion of axial-vector diquarks 
substantially reduces the asymmetries. 
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. , The Sivers function was introduced for the first time in Ref . , in an attempt to explain the observation of single- 
spin asymmetries in hard hadronic reactions. Since then, some phenomenological extractions of the Sivers function 
|, have been performed |^ Q , from data on pion production in proton- proton collisions 0| . The Sivers function 

■ gives a contribution also to the single-spin asymmetry observed by the HERMES collaboration in pion production 
^ [ via deep inelastic scattering off polarized targets 0,I3,I3- In all the above cases, however, the presence of competing 

. effects (in particular the Collins effect) did not allow clear conclusions up to now 0,0. 

' Despite the phenomenological indications, for several years the Sivers function was believed to vanish due to time- 

^ reversal invariance However, this statement was contradicted by an explicit calculation by Brodsky, Hwang 

' and Schmidt, using a spectator model |12|. As the Sivers function is an example of T-odd entity, it requires the 

I interference between two amplitudes with different imaginary parts 0, 0| . Spectator models at tree level cannot 

■ provide these nontrivial phases, but they can arise as soon as a gluon is exchanged between the struck quark and 
_ the target spectator |0. More generally, the presence of the gauge link, which insures the color g au g e invariance 
. of parton distributions, can provide nontrivial phases and thus generate T-odd functions 0, 0, MM The mam 

- ■ ingredient of the model calculation of is nothing else than the one-gluon approximation to the gauge link. It is 



I. INTRODUCTION 



Oh also interesting to note that T-odd distribution functions vanish in a large class of chiral soliton models, where gluonic 
^ ' degrees of freedom are absent (ISJ . 
'T"! [ The work of Brodsky, Hwang and Schmidt was not aimed at producing a phenomenological estimate. A step 



I> ■ forward in this direction has been accomplished in Refs. |19ll20l |. In our article, we present an alternative calculation, 
kjpj ' using the version of the spectator model presented by Jakob, Mulders and Rodrigues in Ref. J2lj. In particular, we 
. include in the model a dynamical axial-vector diquark as a possible spectator, and we explore the Sivers function for 
' down quarks. The necessity to include axial-vector diquarks is also discussed, e. g., in Refs. [23, li^,!!^- Finally, we 
point out that a calculation of the Sivers function in the MIT bag model has been recently presented in Ref. [23 • 

II. UNPOLARIZED DISTRIBUTION FUNCTION /i 

The unpolarizcd distribution function /i can be defined as 

/i(x, 4) = i Tr [{^x, kT;S) + ^x, fcr; S)) 7+] , (1) 
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FIG. 1: Tree-level and one-loop diagrams for the specator-model calculation of the Sivers function. The dashed line indicates 
both the scalar and axial-vector diquarks. 



where S is the spin of the target. The correlator <I>(x, kx) can be written as [T^l 

' , (2) 



^x,K;S) = / e+'^-«(p,5|^(o)/:[o-,oo-]/:[o,,oo,]/:[oo,,«,]/:[oo-,c-]^(OI^,^> 



where the notation C[a,b] indicates a straight gauge hnk running from a to b. In Drell-Yan processes the link runs 
in the opposite direction, to — oo [T^ . For the calculation of the unpolarized function /i the transverse part of the 
gauge link does not play a role and the entire gauge link can be reduced to unity. Therefore, for this first part of the 
calculation it is sufficient to consider only the handbag diagram. 

At tree level, we follow almost exactly the spectator model of Jakob, Mulders and Rodrigues [IJ. In this model, 
the proton (with mass Af) can couple to a constituent quark of mass m and a diquark. The diquark can be both a 
scalar particle, with mass Mg, or an axial- vector particle, with mass A/„. The relevant diagram at tree level (identical 
for the scalar and axial- vector case) is depicted in Fig. ^ (a) . In our model, the nucleon-quark-diquark vertices are 



Ts=9s{k% T^; = ^757^ (3) 



We make use of the dipole form factor 



n2, .r (fc' - (1 - X)2 

9s/vik^) = ^sA, .-.^ '\ .2 : (4) 

+ ^ / J 



where 



4 = -(1 - k^ -xM^/^ +x{l-x) M^, (5) 
Ll/^ ^{l-x)A' + x ~x{l^x) M\ (6) 

The only difference with respect to Ref. j21| is the form of - the vertex involving nucleon, quark, and axial-vector 
diquark. This change modifies the original results only slightly. Note that our choice of the form factor, defined 
in Eq. Q), is very different from the Gaussian form factor employed in Ref. T^. Both choices have the effect of 
eliminating the logarithmic divergences arising from fc^ integration and suppress the influence of the high kx region, 
where anyway perturbative corrections should be taken into account po| . 
The final results for the unpolarized distribution function /i are 

p gl[{xM + mf + k?t\ _ Nni~xf[{xM + mf + kl] 

^^^""'"^^ 2{2^Y{l^x){k^-m^Y 167r3(4 + L2)^ ' 

v( gl[{xM + m)'^ + k'^ + 2xmM] {I - xf [{xM + mf + k"^ + 2xmM] 

^'^""'^^^ ^ 2(2.)3(l-x)(P-m2)2 = 16.3 {kl + LlY ' 
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FIG. 2: Model calculation of xfi{x): with scalar diquarks only (dashed line), with scalar and axial- vector diquarks (solid line). 
The d quark distribution is zero when only scalar diquarks are used. 



Both functions can be integrated over the transverse momentum to give 

(1 - xr 



96 7r2 L6 
96 7r2 L6 



[2{xM + mf + Lf\ , 
[2{xM + mf + Ll + AxmM] . 



(9) 
(10) 



In order to obtain the distribution functions for u and d quarks, we use the following relation, coming from the 
analysis of the proton wave function. 



S 1 

Ji ^ 2 ■' 1 2 



rd j-t) 

Jl — Jl- 



(11) 



Here we refrain from discussing the choice of parameters of the model and its quality, for which we refer to the 
original work ■ We choose the following values for the parameters of the model: 



m = 0.36 GeV, 
A = 0.5 GeV, 



Ms = 0.6 GeV, 
iV^ = 6.525, 



My = 0.8 GeV, 
= 28.716 . 



(12) 
(13) 



The factors Ns and Ny are chosen in order to normalize the functions f( and fi to 1 and consequently to normalize 
/" to 2 and /f to 1. The results of the model are shown in Fig. [5] The dashed line represents the result of the 
spectator model with scalar diquarks only (with /" = 2fl). As can be seen, the difference for the u distribution is 
not big, but it is particularly relevant at small x. The d quark distribution is zero when only scalar diquarks are used, 
which is clearly unrealistic. 

One of the problems when trying to match the model and the phenomenology is that it is not clear at which 
energy scale the model should be applied. A way to estimate this energy scale is to compare the total momentum 
carried by the valence quarks in the model and in some parametrization !24|^3|- Taking, for instance, the GTEQ5L 
parametrization it turns out that this scale is about 0.078 GeV^. Then, by applying standard DGLAP equations, 
we can evolve our model results to 1 GeV^ and compare it with the GTEQ5L parametrization at that scale. The 
result is shown in Fig. 13 Admittedly, the model reproduces the parametrization of the valence quark distribution 
only qualitatively. In any case, in this work we mainly aim at giving rough estimates of the relative magnitude of 
the u and d Sivers function and of the related single-spin asymmetries, as well as studying the changes induced in 
the model results when an axial-vector diquark is introduced. Therefore, we refrain from improving and tuning the 
model. 



^ We use leading order evolution with three flavors and A^i = 0.222, in order to match the CTEQ5 results. 
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(a) (b) 
FIG. 3: Model calculation of xfi{x) (solid line) compared to the CTEQ5L parametrization [g^] (dashed line) at 1 GeV^. 



III. SIVERS FUNCTION AND ITS MOMENTS 



We use the following definition of the Sivers function 

'^^^ fM-, 4) = - J Tr [($(x, fcr; S) - <I>(x, kr; -S)) 7+] , (14) 

where Aiel^kxiSTj = Tr [757+7^ At tree level the Sivers function turns out to vanish. This is due to the lack 
of any final state interaction that can provide the imaginary parts necessary to generate T-odd functions. We need 
to introduce the one-loop amplitude described in Fig.n(b). This is nothing else than the one-gluon approximation 
to the gauge link included in Eq. (j^J. The Sivers function receives a contribution from the interference between 
amplitude (a) and the imaginary part of amplitude (b). The imaginary part of amplitude (b) can be computed by 
applying Cutkosky rules or, equivalently, by taking the imaginary part of the propagator 1/(1'^ + ie) (ISj. Note that in 
Drell-Yan processes the different topology of the one-gluon diagram implies that the imaginary part of the propagator 
1/ — ie) has to be taken, with the effect of changing the overall sign of the Sivers function [ij, [23 ■ This is consistent 
with the change in direction of the gauge link mentioned before |l7| . 

Following Ref. we perform the calculations initially with Abelian gluons and generalize the result to QCD at 
the end. We use Feynman gauge. In order to compute the one-loop diagram, we have to make the appropriate choice 
for the vertex between the gluon and the scalar or axial-vector diquark. We choose the following forms 

= -ie2{2p ~ 2k ~ if , (15) 

r(;'"'3 = [(2p -^2k~ lyg"" - {p - k - vlf g'^'^ ~ {p - k - {1 - i;);) V] , (16) 

where 62 denotes the color charge of the diquark, which we assume to be the same for both kinds of diquark. The 
gluon-axial- vector diquark coupling is identical to the photon-axial- vector diquark coupling suggested in Ref. |22j| . 
The parameter v is the anomalous magnetic moment of the axial-vector diquark; for v — 1 the vertex is analogous, for 
instance, to the standard photon-W vertex. In any case, our results at leading order do not depend on the anomalous 
magnetic moment of the diquark. 

The final results for the Sivers function are 

+ (17) 

^"■'•'•'^'--8(2.).L;[4 + i;l'- ' ' 

Note that our result for the scalar diquark has the opposite sign compared to similar computations 120*1. However, a 
sign error in those computations has been recently identified (291] . The other differences between our results and those 
m Refs. [Hiiiisg are due to the different choice of form factors in the nucleon-quark-diquark vertex. Unfortunately, 
we cannot evolve our results to a higher energy scale, as we have done for the unpolarized distribution function, since 
the evolution equations for the Sivers function have not been established yet ,30„,31,]. We hope, however, that the 
dependence of our results is not very strong. 
We introduce the kj- moments 

f^^'/'\x) ^ I A^kr ^ 4), ^ I d^fc^ _|_ (19) 



5 




(o) (b) 



FIG. 4: Model calculation of xf^^^^\x): with scalar diquarks only (dashed line), with scalar and axial- vector diquarks (solid 
line). The d quark distribution is zero when only scalar diquarks are used. 



for which we get the following results 

6162 (1 - x)^ {xM + m) 



„±(l/2)s 
J IT 



[x) 



fti"'>{x) 



1024 7r2i5 
eie2N'^{l-xfxM 
2048 7r2 Ll ' 



„±(i)i 

J IT 



ix) 



6162 N'^ (1 - xf {xM + m) 
256 7r3ML4 ^ 

6162 jV^^ (1 - xfx 

512n^Lf, 



(20) 
(21) 



The only parameter to be fixed is the product of the quark and diquark charges. Following Ref. [12| we fix 6162 = 
An Cf cts and we choose Cf = 4/3 and as ~ 0.3. 

Relations equivalent to those of Eq. Hll|) hold for the Sivers function and its moments. In Fig. 21 we show the model 
results for the first moment of the Sivers function. The inclusion of the axial- vector diquark results in some change to 
the u Sivers function and allows us to produce a nonzero d quark Sivers function. It turns out in particular that the 
d Sivers function is much smaller than the u one and has the opposite sign. This result is in qualitative agreement 
with the bag-model calculation of Ref. [2^. The opposite sign of the two functions is also in agreement with the 
only phenomenological extractions available at present 0,|3,I13- However, there is a sharp difference in the relative 
magnitude of the two contributions, since in the phenomenological extractions the absolute value of the d contribution 
is about half of the u contribution, while in our model calculation the d contribution is only about 1/10 of the u 
contribution. We point out that this difference could be due to a sizeable contribution of sea quarks (in particular u) 
to the asymmetry studied in Refs. 

Note that our model calculation complies with the positivity bound [131 



(22) 



We performed also the calculation of the function h^, introduced by Boer and Mulders in Ref. |3J| and defined as 



'^j^ hi{x, k^) 

For the scalar diquark, we found that h^^ 
axial-vector diquark, we obtain 



1 



Tr [($(x, kr; S) + $(x, kr; -S)) i a'+Ts] 



(23) 



4 

/j^, confirming the results already obtained in Refs. [T^ |20| . For the 



h-^ (x , krj^ ) — 



6162 N"^ (1 - xf M {2xM + m) 
4(27r)4L2 [kl + Llf ' 



uMi)v( . ^ 6162 (1 - x)^ {2xM + m) 



IV. SINGLE SPIN ASYMMETRIES 



We consider the weighted transverse spin asymmetries 

/ dy d(l>s d^Ph± sm{(j)i, - 0s) (dV[/T - dV[/i) 



(sin (0,1 - (l)s))uT (^'^) 



/ dy d(t>S d^Ph±id'^au^ + d^aui) 
J dy d0sd'A± \Ph±\/izM) sin(0;, - 0s) (dV^/t - dV^/x) 
/d2/d0sd2A±(dVc/T+ dVc/x) 



(25) 
(26) 



FIG. 5: Description of the vectors and angles involved in the Sivers asymmetry measurement. 
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FIG. 6: Model estimate of the single-spin asymmetry {sin {<j)h — 4's))ut- with scalar diquarks only (dashed line), with scalar 
and axial- vector diquarks (solid line). The x and z dependence is shown for tt"*" and 7r~. 



where the notation dtTfji indicates the cross section with an unpolarized lepton beam off a transversely polarized 
target. The angles involved in the definition of the asymmetry are depicted in Fig. [S] These asymmetries are 
currently measured by the HERMES and COMPASS collaborations |3ai36i|. 

Under the assumption that the pion transverse momentum with respect to the virtual photon is entirely due to the 
intrinsic transverse momentum of partons, i.e., Phi_ = zkr, the first asymmetry can be written as 



(sin((/)^ - (j)s))uT (^'^) 



1/^ Eae^/i 



±(l/2)a 
IT 



{x)Dl{z) 



1/^ Y.a^lfl{^)D't{z) 



(27) 



with a indicating the quark flavor. Our model results are displayed in Fig. El We took the unpolarized fragmentation 
functions from Ref. [33 at a scale = 1 GeV^. In order to make predictions useful for the HERMES experiment, to 
produce Fig.|ni(a) and (b) we integrated the asymmetries over z from 0.2 to 0.7. To produce Fig. El(c) and (d) we 
integrated the asymmetries over x from 0.023 to 0.4. 

Evidently, there is not a big difference between 7r+ and ■k~ asymmetries. We do not plot the results for tt*^ 
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FIG. 7: Model estimate of the single-spin asymmetry { ' sin (0^ — (j)s))uT- with scalar diquarks only (dashed line), with 
scalar and axial- vector diquarks (solid line). The x and z dependence is shown for tt^ and 7r~. 



production, since they lie between the previous two. We point out that assuming 



/{'Df"*'") «4/ri^;'^"*'°\ (28) 



the above asymmetry can be written as 



(sm (0. - <^s))uT - 2 + 2 [I + 2 <^ j ' ^''^ 

The first term in this equation is the dominant one. This explains why there are only small differences between the tt"*" 
and 7r~ asymmetry. However, the axial-vector contribution to fi cannot be neglected. From Eq. I|29|l it is also evident 
that the dependence of the asymmetries on z is due to the influence of the unfavoured fragmentation functions. In 
the case of tt" the last term in Eq. H29() is bigger and therefore the z dependence is stronger. In the case of 7r° - not 
shown in our pictures - the term in parentheses would be exactly 3/2 and the z-dependent asymmetry would be a 
flat line at about 5%. 

The asymmetry in Eq. H26|l can be written in an assumption-free way as 

\Ph±\ . .A,^ 1/^ Y.aelfti'^\x)zDt{z) 



The calculated asymmetries are shown in Fig. {7\ where we performed the integrations over x or z as in the previous 
case. As before, assuming Eq. (|28f) the asymmetry can be simplified to 



V. CONCLUSIONS 



We calculated the Sivers function in a spectator model of the nuclcon with scalar and axial- vector diquarks. The 
final state interaction necessary to generate T-odd distribution functions was provided by gluon rescattering between 
the struck quark and the diquark. 
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The inclusion of axial- vector diquarks allowed us to calculate the Sivers function for the d quarks. The function 
turns out to have the opposite sign compared to the u quarks and to be much smaller in size. The u quark Sivers 
function is substantially reduced by the axial- vector contribution. Although the reliability of the model is very limited, 
we think that our results on the relative behavior of u and d quarks could be qualitatively relevant. 

Using the results of our model, we estimated some single spin asymmetries containing the Sivers function. These 
asymmetries arc; at present being measured by the HERMES and COMPASS collaborations. We notic;c(l that the 
inclusion of axial-vector diquarks can make drastic changes in the asymmetries as compared to the spectator model 
with scalar diquarks only, particularly at low x. We were able for the first time to estimate the Sivers single spin 
asymmetry in 7r~ and tt" production. We observed that the tt"*" and 7r~ asymmetries are not very different, due to 
the dominance of the u quark contribution in both cases. The 7r° asymmetries (which we did not show) lie between 
the 7r+ and tt" estimates. 

The present approach does not take into account sea quarks. Unfortunately, at the moment we have no indication 
about the size and sign of the sea-quark Sivers function. 
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